Itinerant electrons in the Coulomb phase 
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We study the interplay between magnetic frustration and itinerant electrons. For example, how does the 
coupling to mobile charges modify the properties of a spin liquid, and does the underlying frustration favor 
insulating or conducting states? Supported by Monte Carlo simulations, our goal is in particular to provide an 
analytical picture of the mechanisms involved. The models under considerations exhibit Coulomb phases in 
two and three dimensions, where the itinerant electrons are coupled to the localized spins via double exchange 
interactions. Because of the Hund coupling, magnetic loops naturally emerge from the Coulomb phase and 
serve as conducting channels for the mobile electrons, leading to doping-dependent rearrangements of the loop 
ensemble in order to minimize the electronic kinetic energy. At low electron density p, the double exchange 
coupling mainly tends to segment the very long loops winding around the system into smaller ones while it 
gradually lifts the extensive degeneracy of the Coulomb phase with increasing p. For higher doping, the results 
are strongly lattice dependent, displaying loop crystals with a given loop length for some specific values of p, 
which can melt into another loop crystal by varying p. Finally, we contrast this to the qualitatively different 
behavior of analogous models on kagome or triangular lattices. 



I. INTRODUCTION 

The combination of magnetism and itinerant electrons is 
a multi-faceted field in the physics of correlated electrons, 
where our understanding is still remarkably patchy: even in 
the case of a square lattice Hubbard model, we lack consensus 
on a detailed phase diagram in the doping-temperature plane. 

Besides the cuprate superconductors, there are plenty of 
other settings in which interesting questions arise, not least 
popularized of late by questions raised by the pnictide super- 
conductors, where magnetic frustration and accidental degen- 
eracies have started to be considered. 

More broadly, there has been increased interest in the in- 
teraction of frustrated magnetism with itinerant electrons^.!!. 
Here, we take up the spirit of this thread of work and study 
itinerant electrons on a highly frustrated lattice, the pyrochlore 
lattice. We consider both three and two dimensions, the lat- 
ter case also being known as the square lattice with crossings, 
planar pyrochlore, or checkerboard lattice. 

We start with an exotic frustrated phase of a magnetic insu- 
lator, the Coulomb phase, which has been extensively studied 
recently^—. This phase has a number of unusual properties, 
including algebraic spin correlations and the emergence of ex- 
tended one-dimensional degrees of freedomi^S, the nature of 
which is an independently interesting problem 2 -^ 

Especially the latter will play an important role in the fol- 
lowing analysis, given the ability of electrons to provide ev- 
idence for non-local structures through properties related to 
transport phenomena. Indeed, it is the marriage of the lo- 
cal constraints imposed by frustration with the 'non-local' 
physics describing mobile particles which makes up for much 
of the interest in this field. 

The Ising-double exchange model which we study here has 
many parameters: electron density, p, temperature T, Ising 
anisotropy, Hund's coupling, Jh, magnetic exchange, J, and 
electron hopping integral, t. A full study of general parameter 
choices is well-nigh impossible in any detail analytically. 



Our approach to the problem considers a regime where the 
effects of frustration are particularly strong but where consid- 
erable progress towards a detailed description is nonetheless 
possible by analytical (or simple numerical) means. That a 
(non-trivial) regime where this is possible exists at all is a pri- 
ori not obvious, and we find that we need to restrict a number 
of parameters to limiting values: we study the limit where the 
magnetic energy scales are much larger than the hopping inte- 
gral t, so that the resulting problem is one of electrons hopping 
on a classical background spin configuration. We do not need 
to restrict the electron density to be small, although for that 
case, we have the most detailed set of results. 

We find that the low-density behaviour can be mapped onto 
a study of a classical loop model with non-trivial weights aris- 
ing from the addition of electrons Hund's coupled to the spins 
in the Coulomb phase. This results in phenomena such as a 
transition from a (in two dimensions, critical) percolation sit- 
uation to one in which the loops acquire an exponential length 
distribution, thereby removing all conducting paths across the 
sample. As the doping is increased further, we find a sequence 
of density-dependent preferred loop lengths, which lead to a 
tendency to form loop crystals which may, however, be frus- 
trated by the lattice geometry. 

The organization of this article, and our main results, are 
summarized below. 



A. Summary and Overview 

We restrict ourselves to the limit of large exchange cou- 
pling interactions so that the ice-rules (see figure [TJ them- 
selves are never compromised. We make extensive use of the 
loop picture encoding the ice rules.-- The loops serve as ID 
channels for the electrons. The problem is thus transformed 
to entropy and energy considerations of possible loop cover- 
ings, with loops supporting varying numbers of electron. This 
"loop framework" for describing conduction electrons is in- 
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traduced in SectionHIl 

Ref. llO has described in some detail the loop distributions 
for both 2D and 3D in the absence of electrons. Once electrons 
are added, within each loop electrons can occupy states whose 
energy is given by a ID dispersion. The dispersion minimum 
is independent of loop length, so the first electron in a loop 
has the same energy in all loops. As a result, there is a low 
doping regime where it is possible to fit at most one electron 
per loop. All such configurations have the same energy. This 
is the entropic regime, because entropic arguments determine 
favorable configurations within an equal-energy manifold. 

When the density of electrons p is larger, a sub-extensive 
number of states, possibly even a unique one, tend to be fa- 
vored, because they manage to minimize the kinetic energy; 
this we term the energetic regime. 

In Section [TTTJ we present entropic considerations relevant 
to the regime of low doping. The total entropy contribution 
comes from both the loops (with parameters extracted numer- 
ically), and the electrons (derived analytically). In particular, 
the presence of electrons acts like a cutoff on the total number 
of loops in the system. In 3D, forbidding configurations with 
less loops than electrons suppresses the formation of exten- 
sive loops (present in the peak of the probability distribution 
function (PDF), see e.g. figure [5]), but does not modify the 
exponents of the PDF. On the other hand in 2D, there is evi- 
dence for a variation in the exponent of the power law of the 
loop distribution. 

Once the electronic density gets too large for loops to be re- 
stricted to at most one electron, we need to consider energet- 
ics. Section HVl presents these energy considerations and the 
phase diagram obtained thereby. Ignoring lattice constraints 
on loop coverings, we use energy calculations and a Maxwell 
construction to obtain the phase diagram of Figure[9] We ver- 
ify some of these results through Monte Carlo simulations. 
We also present constraints imposed by the lattices under con- 
siderations and identify loop crystals arising as the doping is 
varied. 

With periodic boundary conditions, loops spanning the sys- 
tem can be divided into segments connecting "opposite" faces, 
which we call filaments. Because the transmission of elec- 
trons through the system can only occur via these conduct- 
ing channels, section [V] is dedicated to their statistics, as a 
function of doping p and dimension. In absence of itinerant 
electrons, the number of filaments grows linearly with (cubic) 
system size in 3D but remains constant and of 0(1) in 2D. 
While these behaviors are qualitatively not modified at low 
doping, the conducting channels turn out to vanish at interme- 
diate values of p. 

A separate final section is devoted to an outlook which also 
contains some words on the behaviour of analogous models on 
the triangular and kagome lattices, which turn out to exhibit 
qualitatively different properties both from the pyrochlores 
and from each other: we find a magnetic conducting solid as 
well as an insulating cooperative paramagnet. 

The numerical component of our work involves Monte 
Carlo simulations of several types; some details are provided 
in the Appendix. 




FIG. 1: In the Coulomb phase, each frustrated unit (crossed squares 
in 2D checkerboard, left, and tetrahedra in 3D pyrochlore, right) pos- 
sesses two up and two down spins, respectively colored in blue and 
red. Those are the so-called ice-rules or divergence free conditions. 
Connecting spins of the same color forms a network of loops, as illus- 
trated for the checkerboard. This model is equivalent to the nearest 
neighbour spin ice model. 



n. THE SYSTEM AND THE LOOP FRAMEWORK 
A. The model 

We will focus primarily on the checkerboard and py- 
rochlore lattices, which are two lattices where a Coulomb 
phase can appear (see figure [TJ. The localized magnetic mo- 
ments are Ising spins Si all parallel to a global axis, while 
itinerant electrons can hop on the lattice sites. The Hamilto- 
nian is 

(»>i> 

- Jh ^ 4,a( a a,P ■ Si) djj (1) 

where t is the hopping integral between two neighboring sites, 
c\ a (ci, a ) are creation (annihilation) operators of itinerant 
electrons of spin a on site i, and a a ^ are the Pauli matrices. In 
order of appearance, the terms in equation ([TJ are the antifer- 
romagnetic nearest-neighbour exchange between the localized 
spins incorporating magnetic frustration, the hopping term al- 
lowing movement of itinerant electrons, whose spins interact 
with the localized magnetic moments through ferromagnetic 
Hund coupling (last term). 

In this work, we focus on the limit t <C Jh "C J ■ In 
this limit, the highly degenerate ground state of the frustrated 
system serves as background for the motion of the electrons. 
Magnetic excitations (violations of ice rules) are not present in 
this limit. Electrons can only hop between nearest neighbour 
spins having the same orientation. 

At zero temperature, a Neel or ferromagnetic order would 
give rise to an insulating or metallic state respectively, but a 
spin liquid provides a network of conducting paths for itiner- 
ant electrons. The present work unveils the geometry of this 
network, as it is influenced by both the lattice and the mini- 
mization of the hopping energy. 



B. Emergence of loops 

The pyrochlore and checkerboard lattices are made of 
corner-sharing units with four spins; respectively the tetrahe- 
dron and the square with crossings. The antiferromagnetic 
couplings impose the so-called "ice-rules" with zero mag- 
netization per unit, obtained with two spins pointing up and 
two spins pointing dowr*22. In absence of itinerant electrons, 
this ground state is highly degenerate with 6 possible con- 
figurations per unit. It corresponds to the 6-vertex model in 
2 dimensions 23 and can be mapped onto the nearest neigh- 
bor spin ice model in 3 dimensions 2 *^. These spin systems 
serve as background for the emergent physics of the so-called 
Coulomb phased, a gauge theory where the discrete ice-rules 
under coarse-grained lead to the emergence of a divergence 
free flux. 

Joining spins of the same orientation in every unit, one ob- 
tains loops of up spins and loops of down spins (see Fig. Q]). 
The Coulomb phase can thus be described as an ensemble of 
possible loop coverings. The resulting loop model possesses 
two flavors (loops of up and down spins), where every site of 
the premedial lattice^ is occupied by two loops, one of each 
flavor, and every bond is visited by one loop only. In previous 
work by some of the authors^, a detailed account has been 
given of the statistics and distribution of these loops, both for 
the 2D and 3D cases. 

In the limit of large J/ t , all up (down) electrons are only 
allowed to hop along an up (down) loop, and are constrained 
to remain within this loop. This reduces the electron dynamics 
to be one-dimensional whatever the dimension of the lattice. 

The ID hopping restriction allows us to describe electron 
dynamics in terms of the dispersion of a ID tight-binding 
problem, Ek = — 2tcosk, with k the ID momentum along 
the loop. The ID momentum along the loop, k = 2nq/l, is 
discrete for a loop of finite length I (q = —1/2, 1/2 — 1). 
An up (down) loop of length I can contain between to i up 
(down) electrons. Note that double occupancy does not occur 
in the limit we are considering. 

The lowest and highest single-particle levels in the disper- 
sion have energy ±2t, independent of the loop length (fig- 
ure 

Our system is thus described by 

• the number of lattice sites TV, equal to 4 L 2 and 16 L 3 
for the checkerboard and pyrochlore lattices respec- 
tively, where L is the linear number of unit cells; 

• the total electron number N e , and the electron density 

P = N e /N; 

• the loop histogram of a given configuration, i.e. the 
number hi of loops of length ti, ii being necessarily 
an even number on a bipartite lattice; 

• the type of lattice, which will among other things de- 
termine the smallest possible loop in the system l m in 
(4 for checkerboard and 6 for pyrochlore); the longest 
possible loop length is always ~ N/2. 

We define a few additional relevant observables: 
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FIG. 2: Single-particle energy levels of itinerant electrons confined 
to loops of length I = 8 (left) and £ = 12 (right), due to the hopping 
term in Hamiltonian l[T). There are nondegenerate levels at the high- 
est and lowest energies, ±2t, independent of loop length; the other 
levels are doubly degenerate. 

• the total number of loops in a given configuration Ne = 

• the average loop length for a given configuration 

j = Si Mi = N_. 

• the statistical average loop length (£) over all loop con- 
figurations; 

• the statistical average number of loops (Ne); 

• the number of filaments (see section Ml. 



III. LOW DOPING REGIME 

In this section we consider the low doping regime where N e 
is small enough to have loop configurations with more loops 
than electrons, Ne > N e . Since the lowest single-electron en- 
ergy level in any loop is —2t, the minimum accessible energy 
is the same (—2tN e ) for all such configurations. Therefore, 
the ground state manifold consists of all such loop coverings 
with the same energy —2tN e . The free energy within this 
manifold is then only determined by entropies. At zero tem- 
perature, entropy is understood in the sense that all configura- 
tions have an equal probability to occur. 

In the first subsection below, we present some entropy cal- 
culations, combining loop and electronic contributions to the 
entropy, and show how this determines the average loop length 
at nonzero electronic density p. In the second subsection we 
present numerical results on the effect of electrons on the en- 
tire loop length distribution (loop PDF). The effect on the loop 
PDF is a natural way to characterize the influence of electrons 
in the magnetic Coulomb-phase system. 
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FIG. 3: Histogram of distribution of the number of loops Nc per configuration in 2D (left) and 3D (right), in the absence of electrons. The 
distribution is scaled according to the Gaussian expression of equation ([2} collapsing all system sizes onto the same curve (L — 100 to 600 in 
2D and L — 4 to 60 in 3D). In 2D, the two Gaussians have the same width and corresponds to even and odd Ne (upper and lower curve). In 
3D, finite size corrections are visible for L = 4 (red crosses). 



A. Entropy 



loop contribution to the entropy (~ In P) is 



The total entropy consists of loop and electronic contri- 
butions. The loop distribution is of course itself affected by 
the itinerant electrons. However, in the limit of small p, the 
change in loop entropy is small. Below we combine the p = 
loop entropy with the finite-p electronic entropy to approxi- 
mate the total entropy at small p. 

Loop entropy. Figure 3 shows the distribution of the num- 
ber of loops Ni having length I, in the absence of electrons 
(p = 0), obtained from Monte Carlo simulations (Appendix 
lAl) . The distribution has Gaussian form: 

(N e -(N e ))^ 



P(N e ) 



1 



exp 



2kN 



(2) 



where n 2 d = 0.0384 and K 3d = 0.00423. In 2D, there are two 
Gaussians, corresponding to even and odd Ni. 

For a system with short-distance correlations, Gaussian dis- 
tributions are natural to expect from the Central Limit Theo- 
rem, since by dividing the system into small mesoscopic seg- 
ments the total distribution can be recast into a sum of many 
random variables. In our case, however, we have a system 
with algebraic correlations and extended objects (loops), so 
finding Gaussian distributions is not a priori trivial. 

Equation [2] can be expressed in terms of the average loop 
length I = N/N e instead of Nf. 



P(£) = P(N e ) 



df 



exp 



TV 
2k 



exp 



N 
2k 



1 

1 ~ 



(3) 



Here we have used !/£{£) ~ l/(£) 2 , which is valid for large 
N in the region where the Gaussian is appreciable. Thus the 



Sloop — Si 



2k (I) 



(4) 



where S\ is a constant. 

Electronic entropy. Since we have at most one electron 
per loop, the number of possible combinations to put N e 
electrons in Ni loops is the binomial (Pascal) coefficient 
Ni\/[N e l(Ni — N e )\]. The logarithm then gives the electronic 
contribution to the entropy. Using Stirling's approximation 
for the thermodynamic limit (Ni >• 1, N e ^> 1), we get as 
per usual 



N 

Seiec = j[x\n.x + (l-x) ln(l - a:)] 



(5) 



where x = N e /N l = pi. 

Total Entropy. For small electron densities, we can ap- 
proximate the loop entropy by the p = expression calculated 
above. The total entropy of electrons and loops is then 

Stot = -j[p£]npl+(l- pi) m(l - pi)] 

N - ? 
+ S 1 -^ m (l-(t) ) . (6) 

The subscript represents the ensemble average at p = 0. The 
value of 1 that maximizes Stot is found (in the limit pi -C 1) 
to be 



io P t w (£) - P k{£)1 . 



(7) 



For p(£)o <C 1, the distribution stays almost Gaussian, so this 
most probable value of £ is approximately the mean value of 
the distribution (I) p . In Figure [4] we compare this prediction 
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FIG. 4: Mean value of the loop length (I) as a function of doping. Left: checkerboard, L — 10 (red dots) and 60 (blue triangles). Right: 
pyrochlore, L = 6 (red dots) and 20 (blue triangles). The dashed straight lines are predictions from Equation Q, without any fitting parameters. 
For both panels, the blue triangles data points are very close to thermodynamic limit. 



to numerical data. For each system size L, we can extract (£)q 
from the numerical (£} at p = 0. Eq. © then gives a linear 
prediction (dashed straight lines in Figure [4}, which works 
well for small p. 



B. Loop length distribution 

We next examine the effect of electrons on the entire loop 
probability distribution function (PDF). We denote by r the 
power-law exponent, when the PDF has form £~ T . Without 
electrons (p = 0), the PDF follows P 2D ~L 2 /£ T in 2D with 
t = 2 + 1/7. In 3D, the p = PDF displays a crossover 
around I ~ L 2 between two power laws, from L 3 /£ 5 / 2 to 
1/1 (For details, see Ref.|H) 

At finite electron densities, loop configurations with Nt < 
N e are rejected due to energetics, as explained previously. 
Thus electron doping acts not entirely unlike a "chemical po- 
tential" for loops, favoring configurations with more loops, 
and thus a priori shorter ones. In both 2D and 3D, this implies 
a disappearance of longer loops, as can be seen in the calcu- 
lated PDF's of Figures [5] and [6] where the large-^? parts of the 
PDF are progressively decimated for increasing p. In 3D the 
form of the PDF is otherwise unchanged (Figure|5]l. In 2D the 
effect seems to be more drastic; the inset of figure [6] suggests 
that the entire power-law behavior of the PDF is modified. 

To quantify how the 2D loop PDF changes qualitatively 
at finite p, we define and compute a local exponent in I: 
riocai(^p) = \og{P 2D (i,p)/P2D(2£,p)) /log2. This is dis- 
played in the upper right inset to Figure[6] Our results suggest 
a trend toward increasing ri oca i as a function of p, consistent 
with the idea that itinerant electrons favor small loops. 

This outcome deserves a few comments. In Ref. |2(| the 
loop statistics of the 2D Coulomb phase (zero doping) has 
been shown to be analogous to the Stochastic Loewner evolu- 
tion process SLE K= 6 with fractal dimension Dj = 1 + re/8 = 
7/4. The SLE K can be identified to various realizations of 
the 0(n) model through the relation n = — 2cos(47r/re)i2& 
The 0(n) model is often used to describe fully packed loop 
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FIG. 5: Loop length distribution Pzd{£, p) in 3D for different val- 
ues of doping (p = {0, 0.00075, 0.0015, 0.002}]) and system size 
L = 10. The distribution is normalized such that J IPzoil, p)A£ = 
16L 3 = N. The dashed line indicates the power law fit £~ 1 at p — 
and the arrow shows the shifting of the distribution for increasing 
p. The exponents of the two power law regions (before and after 
I ~ L 2 = 100) do not vary, but the peak for long winding loops 
t ~ 8L 3 = N/2 gets smaller with increasing p. 



models with loop fugacity n (see e.g. Ref. |23|) . The partition 
function of the fully packed loop model is Z = ^ n Ne , where 
the sum runs over all possible configurations. The Coulomb 
phase corresponds to a fully packed loop model™; at zero dop- 
ing, the free energy of our model is trivially independent of the 
number of loops Nt and thus corresponds to a fugacity n = 1. 
Higher values of the fugacities favor configurations with more 
loops and tend to increase the value of r— , in a way reminis- 
cent of the influence of doping here. The addition of itinerant 
electrons remains a non-trivial problem and is not exactly the 
same as a fugacity for a loop, but at small and finite doping p, 
some features could be captured by 0(n{p) > 1) models or 
SLE K ( p ) <6 processes. 
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FIG. 6: Main: Loop length distribution Pid (£, p) in 2D for dif- 
ferent values of doping (p = {0,0.0025,0.0056,0.0087}]) and 
system size L = 40. The distribution is normalized such that 
J £P-2d{£, p)d£ = 4L 2 — N. The dashed line indicates power law 
L 2 /£ 15/7 and the arrow shows the shifting of the distribution for in- 
creasing p. On this scale, it is not obvious whether the exponent of 
the power law varies or not. Inset: The power-law exponent defined 
locally on I, n oca i for £ = 8(B), 16(a), 32(«) for L = 40 (similar 
behavior obtained for L = 20 and 60). 



IV. LARGE DENSITIES; PHASE DIAGRAMS 

In this section, we present results relevant for higher den- 
sities, where entropic considerations are no longer sufficient, 
and non-trivial electronic hopping energies need to be con- 
sidered. We continue to describe the system in terms of loop 
coverings. 

We first provide an analysis based on calculations for loop 
coverings of equal-length loops. Using a Maxwell construc- 
tion, we can use this information to predict ranges of electron 
density where the ground state manifold consists of coverings 
by loops of two different lengths. 

These considerations, described in the first two subsections 
below, do not take into account any lattice constraints other 
than the fact that the minimum loop length is £ m i n = 4(6) 
for the checkerboard (pyrochlore). Lattice constraints, disal- 
lowing some loop coverings, are difficult to enumerate or list 
comprehensively on account of their non-local nature. In the 
final subsection, we present the lattice constraints that we have 
identified, and their implications. 



A. Equal-length loop configurations 

In this subsection and the next, we imagine that coverings 
with any unique loop length £ ^ £ m in are possible. Later, we 
will show that at certain fillings, the ground states exhibit a 
unique loop length, while at others, the behaviour can be un- 
derstood through a Maxwell construction based on the single- 
length results. 

For an odd number of electrons p£ = 2n Q + 1 in a loop of 



length £, the energy is (see figure |2|i 

„. „. ( 27m \ sin (tto) 

E M = -2t £ cos(— J =-2*^ (8) 

n=—n ^ / \ t / 

For an even number of electrons 2n G , this expression becomes 

E(p, £) = -2t SilW( ;~ 1 / g) - 2t cosM (9) 
sin (7r / 1) 

More generally, if pt is not an integer, we define the highest 
odd integer below p£ as 



(10) 



where E(.) is the floor function. Each loop is filled with at 
least ?/ electrons up to the energy levels at k = ±ir(ri — l)/£, 
while the remaining N(p — r)/l) electrons in the system are 
distributed in the partially filled level at k = ±tt(ii + l)/£. 
The total energy is 



E(p,£) = -2tN 



sin (irr]/£) 
£sm(n/£) 



1 



f 



11) 



These expressions are electron-hole symmetric, i.e. invariant 
under p f-> (1 — p). 

Figure [7] displays the loop length l(p) that minimizes the 
energy ( fTTT i as a function of p, and thus corresponds to the 
ground state if we impose a unique loop length in the sys- 
tem. (The electron-hole symmetry shows up clearly through 
the mirror symmetry on either side of p = 0.5. Therefore, we 
shall from now on only consider densities below 1/2.) 

Since the lowest possible energy for an electron is — 2t and 
is only accessible for one electron per loop, loops of length 
4 are favored up to p = 1 /4. Of course for lower densities 
than 1/4, other configurations may be possible as long as the 
number of loops is larger than the number of electrons, but a 
system with only loops of length 4 will always be part of the 
ground state manifold for p ^ 1/4. 

Remarkably, for most densities, a finite loop length with 
discrete energy levels is preferred compared to infinite loops, 
except for p = 1/3 and 2/5. 

We plot on figure [8] the minimum energy corresponding to 
the loop length £(p) of figure|7] E mm (p) = E(p, l(p)). 



B. Maxwell construction; lattice-independent phase diagram 

We now move beyond configurations with unique loop 
length. We need to consider mixtures of electron densities in 
the different loops. This is done through a Maxwell construc- 
tion, similar to the physics of a liquid-gas first-order phase 
transition. 

A system of N sites and density p can be divided into two 
subsets of sites Ni and 7V2, with density p\ and P2 of electrons 
and loop length £\ and £2 respectively, with 



N p 



= Ni 
Ni pi 



N 2 

N 2 P2- 



(12) 
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FIG. 7: Loop length £(p) minimizing the energy il lb for a system filled with loops of length I only, assuming no lattice constraints other 
than Imin = 4 (left panel, checkerboard) or 6 (right panel, pyrochlore). p takes all rational values p/q, with q £ {0, 1, 2, 1000} and 
p 6 {0, 1, 2, q}, and we consider all loop lengths from £ m m to £„ la x = 100 (red crosses) and 1000 (blue squares): plotting two different 
values of £ ma x provides a graphical way to visualize the two values of p where the most favorable loop length is infinite, namely p — 1/3 and 
2/5 (and their symmetric images with respect to p — 1/2). If several loop lengths give the same energy, we plot the smallest one. 



Now if a straight line between E m i n (pi) and E m i n {p2) re- 
mains below the curve E m i n {p) on figure [8] i.e. if 

N E mm (p) > NxEmi^px) + N 2 E min (p 2 ), (13) 

then the mixture of two densities is more stable than a unique 
density and "phase separation" occurs. We thus construct 
a phase diagram, separating regions of different loop length 
combinations. 

On figure[9] the dots correspond to densities where a unique 
loop length is favored, whereas the zones between them are 
"phase mixed", defined by the loop lengths and electron 
densities of the surrounding dots. The ratio of one phase 
compared to the other is given by equations (fTZt . For example 
for p = 0.4, 80% of the sites belong to loops of length i\ = 8 
with electron density pi = 3/8, while the remaining 20% 
belong to loops of length £2 = 6 with density P2 = 1/2. 
As shown on figure [7] loops of length 4 are particularly 
robust over a wide range of p; preventing their formation, 
e.g. in the pyrochlore lattice in d = 3 which permits loops 
of minimal length 6, thus strongly modifies the phase diagram. 

These results indicate that itinerant electrons tend to favor 
relatively small loops and to prevent the formation of infinite 
ones at zero temperature. 

At low electron densities (up to l/£ m i n ), the minimum en- 
ergy for a given p is degenerate: e.g. for p = 1/12, having 
only loops of length 4, 6, 8, 10 or 12 gives the same energy, as 
the number of electrons is smaller than the number of loops on 
the system, and every electron can fill the lowest energy level. 
This explains why E m i n (p) is a straight line in this region 
(Figure |8): as p decreases from l/£ m in to 0, the degeneracy 
of the ground state increases until one recovers naturally the 
full degeneracy for p = 0. This of course corresponds to the 
entropic regime mentioned above. 

We stress again that these results have been obtained in- 
dependently of the lattice (except for the value of £ m in), and 



that lattice constraints (next subsection) will modify some of 
the phase diagram. 



C. Lattice constraints on loop coverings 

We now consider effects of the lattice in disallowing some 
of the configurations predicted by our analysis above. For the 
specific densities where it is possible to cover the lattice with 
loops of a unique length, the result is an ordered loop crystal. 
In some cases, however, a covering by a unique length or by a 
combination of loops of two lengths, is not possible. We point 
out some such cases below. 



1. Checkerboard 

We first focus on the regime p € [1/4; 3/8] favoring loops 
of length 4 or 8 according to the Maxwell construction. As 
illustrated in figure [10] not only can the lattice be covered by 
these loops, but the transformation from 4 loops of length 4 to 
2 of length 8 is also purely local and allows all possible ratios 
between these two phases in the thermodynamic limit. A pair 
of loops of length 8 (one made of up spins, the other of down 
ones) cannot be separated if there are no other lengths than 
4 and 8 in the system. We shall call such pair a defect. The 
defect concentration is determined by p. Defects are not topo- 
logical in the sense that they can be created and annihilated 
locally, they can be placed anywhere in the background made 
of loops of length 4. 

At quarter filling, there is one electron per loop of length 
4, its energy is — 2t and their total number is N/4. Each 
additional pair of electrons fills the first excited level of the 
newly created defect (3 electrons per I = 8 loop) and gives 
an energy -At(V2-l). At a given p, the number of elec- 
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FIG. 8: Minimum energy E(p) corresponding to the value of £(p) plotted on figure [7] for a system filled with loops of length £ only, 
assuming no lattice constraints other than l m in = 4 (left, checkerboard) or 6 (right, pyrochlore). p takes all rational values p/q, with 
q e {0, f , 2, 1000} and p G {0, 1, 2, q}. We consider all loop length from £ min to £ max = 1000. We chose arbitrarily t — 1/2. A 
Maxwell construction can be visualized from this figure. The red dots are a set of points, such that the dashed lines connected them are always 
below the red curve E(p). Hence a mixture of two phases corresponding to two consecutive red dots (£i, pi) and (£2, p-i) has a lower energy 
than a single phase with a unique density of electrons p on a unique type of loop length £. For increasing values of p, the red dots correspond 
to the most favored loop length £ = {4, 8, 6, 8, 4} (left panel) and £ = {6, 12, 10, 8, 6, 8, 10, 12, 6} (right panel), as can be read from figure|7] 



4 8 6 I 
Checkerboard • • • • * 

Q 1/4 3/8 1/2 P 

6 12 10 8 ti i 
Pyrochlore • m m • • • ► 

Q 1/6 1/4 3/10 3/8 1/2 P 



FIG. 9: Zero temperature "phase diagram" of loop length configuration as a function of electron density p, obtained without accounting for 
lattice constraints other than £ min =4(6) for checkerboard (pyrochlore). A phase is defined by the length of its constituting loops and the density 
of electrons on them. The values of p and £(p) for each phase are given below and above the lines. 
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FIG. 10: Schematic representation of the checkerboard lattice with 
16 sites and different loop coverings. Left: 4 loops of length 4, cor- 
responding to the ground state at p — 1/4 with 1 electron per loop. 
Center: 2 loops of length 8, corresponding to the ground state at 
p = 3/8 with 3 electrons per loop. Between 1/4 and 3/8, a mixture 
of these two configurations will occur, their respective ratio being set 
by the total number of electrons N e = pN. This arrangement of 
two loops of length 8 appears as a "defect" in a crystal of loops of 
length 4. Right: a red loop of length 6 spanned by another blue loop; 
the shortest way to close this blue loop requires 10 sites. It is thus 
impossible to have a mixture of loops of length 6 and 8 only. 



trons added with respect to the loop crystal at quarter filling is 
{pN - N/4) = (4p - l)L 2 . The total energy of the system 



between p = 1/4 and 3/8 is then 



E(p) = -2tL 2 (l + (4p-l)(V2-l)) (14) 



Once normalized per number of sites 4L 2 , this expression is 
the dashed line plotted on figure [8] We have verified this ana- 
lytical result with finite-temperature Monte Carlo simulations; 
some details of the method are in Appendix [A] Extensive de- 
generacy is recovered in this region between 1/4 and 3/8. 

At higher doping, above p = 3/8, the Maxwell argument 
predicts a mixture of loops of length 8 and 6. However, as 
we can see in figure [10] a single loop of length 6 imposes 
the presence of loops of length 10 at least. This implies that 
this region of the phase diagram (obtained without accounting 
for such lattice constraints) is further modified, in an as yet 
unknown manner. 
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2. Pyrochlore 



Analogous modification of the phase diagram of figure[9]is 
more severe for the pyrochlore, and shows up already at p = 
^-/f-min = 1/6. We have found that it is impossible to cover 
the pyrochlore lattice with loops of length 6 only. However, it 
is possible to do it for half of the system, as explained next. 

The pyrochlore lattice can be seen as an alternative stack of 
kagome and triangular layers orthogonal to one of the global 
[111] axes. As depicted on figure [Di d). 2/3 of each kagome 
layer can be filled with loops of length 6 (blue hexagons), 
while the other 1/3 of the kagome sites forms extensive 
winding loops along the [1 1 l]-axis, crossing alternatively the 
kagome and triangular layers. Since all blue sites form loops 
of length 6, putting one electron per loop provides a ground 
state configuration up to p = 1/12 at least. This is a priori not 
the only one, but this proves its existence. 

Indeed, despite an intensive search by complete enumer- 
ation of configurations respecting the ice-rules on the py- 
rochlore lattice, we have not detected any relevant loop crystal 
or mixture of them for systems up to 128 sites. A system of 
size 192 has been partially investigated, with the same out- 
come. The smallest occurrence of a single-length loop cover- 
ing is for 8 loops of length 16 in a system of 128 sites; however 
this is not a relevant length according to our phase diagram of 
Figure [9] Thus, the loop coverings on an actual pyrochlore 
are modified from Figure|9]in most or all of the density range 
p G [1/6,1/2]. The details of this modification remains an 
important open problem. 



V. CONDUCTING CHANNELS OR "FILAMENTS" 

In this section, we discuss the effect of the Coulomb phase 
and loop structure on the conductivity of mobile electrons. 
Since the electrons are confined to loops, they can conduct 
only if a loop connects one edge of the sample to the opposite 
end. Therefore, we study the number of such sample-spanning 
loops, which, following Ref.|28l we refer to as "filaments". 

When there are no filaments {e.g. in loop crystals with only 
finite-length loops), the system is unambiguously an insula- 
tor. When there are filaments spanning the system, the system 
cannot be immediately called a conductor, because the actual 
conductivity will depend on scattering mechanisms exterior to 
our model. 

In the first subsection below, we consider small dopings, 
where we present Monte Carlo results for the average num- 
ber of filaments as a function of system size. In the second 
subsection, we comment on the consequences of our phase 
diagram. 



A. Filaments at low densities 

We first consider very low densities that the loop distribu- 
tion can be assumed to be largely unchanged from the p = 
case. 
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FIG. 11: Average number of filaments, i.e. those segments of loops 
spanning the entire system from one border to the opposite one, as 
a function of system size L. We use two definitions of filaments; 
they can either cross the orthogonal borders thanks to the periodic 
boundary conditions (•) or not (■). The cartoon on the right shows 
such a boundary-crossing filament, which would be excluded in the 
■ data. Dashed lines are guides to the eye for the linear behavior 
with L. Both x- and y-axes are on a logarithmic scale. 



In 2D, there is a small but constant number (w 1.86) 
of winding loops on average in the thermodynamic limit 20 . 
There is thus some probability to have filaments in the 
checkerboard case, but the number of filaments does not grow 
with system size. 

In 3D, the background of winding loops ensures that there 
are filaments whose number grows with system size. The data 
in FigureQT]shows that the number of filaments increases lin- 
early with the linear size of the cubic sample. 

Since our simulations and loop counting are performed with 
periodic boundary conditions, some of the filaments, while 
spanning the sample in one direction, also cross one of the 
orthogonal boundaries. (An example is shown in the cartoon 
to the right of Figure QT|) One can argue that this type of 
filament would not contribute to conduction in a real-life cubic 
pyrochlore sample. Therefore, we show data both excluding 
and including this type of filament, and they are seen to have 
the same power-law behavior. 

Since the data including these filaments have better statis- 
tics, for finite p we display the inclusive data, with the expec- 
tation that there is no qualitative difference. 

According to Fig.|5j itinerant electrons tend to make disap- 
pear extensive loops of length L 3 in 3D. Whether or not this 
prevents the formation of filaments is not as straightforward as 
it seems. Indeed, one could naively assume that the number 
of conducting channels will decrease and maybe even vanish 
in the thermodynamic limit. However the number of filaments 
remains approximately constant as plotted on figure [12] This 
means that the average number of loops increases with dop- 
ing via dividing the very long ones into smaller but nonethe- 
less extensive loops spanning the system. In the low doping 
regime in 3D, the number of conducting channels is approxi- 
mately independent of p. 

By contrast, in 2D, with relatively large error-bars, the num- 
ber of filaments decreases but remains of O ( 1 ) as electrons are 
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added. 



B. Conduction channels at larger fillings 

In Section llVl we identified densities (for both 2D and 3D) 
where the system is a loop crystal or is covered by loops of two 
finite lengths only. In such cases, there are no filaments, and 
the system is truly insulating: lifting the frustration-induced 
degeneracy removes the non-locality of the loops. How the 
energetics (and resulting degeneracies) imposed by lattice 
constraints (see section IIV Cl > manifest themselves in trans- 
port properties is an intriguing separate question. 

For some ranges of p, especially in 3D, the loop cover- 
ings indicated by the Maxwell construction are disallowed by 
the lattice geometry (section llVCb . In such cases it remains 
an open question whether or not the lattice constraints result 
in loop coverings including infinite (sample-spanning) loops. 
Unfortunately, the issue of conducting channels depends on 
the answer to this generally unresolved question. 

VL SUMMARY AND OUTLOOK 

We have analyzed the double exchange model on the py- 
rochlore lattice in two and three dimensions. We have chosen 
to consider a parameter range for which (i) magnetic frustra- 
tion is known to give rise to unconventional ground state en- 
sembles and (ii) where it is possible to make considerable an- 
alytical progress by mapping the system onto an ensemble of 
loops, the statistical properties of which are influenced by the 
addition of the electrons. 

We have identified a number of phenomena which depend 
on features such as dimensionality, which determines whether 
or not there exist loop segments winding around the system; 
or lattice structure, which may frustrate the geometric packing 
of preferred loop lengths. 

The model studied here leaves unanswered a number of 
questions and immediately suggests many generalizations and 
extensions. We have worked in a limit of parameters such 
that the magnetic exchange and anisotropy dominate over the 
Hund coupling which in turn dominates the hopping integral. 
Our analysis applies to zero temperature. 

It would be interesting to relax any of these choices, al- 
though technically this may not be easy. In particular, given 
the presence of gapless excitations on long loops, interesting 
low-temperature physics may appear. Canting can give rise 
to non-trivial Berry phase physics, and finite Hund's coupling 
will enable electrons to hop between loops. For example, re- 
cent work at quarter filling on the pyrochlore Kondo lattice 
has shown the emergence of a chiral magnetic order in the 
weak-coupling regime^. 

Finally, even in the parameter range discussed here, it will 
be interesting to ask how different frustrated lattices shape 
up compared to the pyrochlores. We devote the final para- 
graphs of this paper to discussing this question for the case 
of triangle-based lattices, the triangular Bravais lattice and the 
Archimedean kagome lattice. 



A. Other lattices: triangular and kagome 

The Ising ground states observed for the pyrochlore lat- 
tice have vanishing total spin on each tetrahedron. These 'ice 
rules' ensure the existence of the Coulomb phase and states 
obeying them amount to a moderate zero-point entropy of less 
than a third of that of a free spin. By contrast, the zero-point 
entropy of the triangular Ising magnet is not far from half of 
that of a free spin, while that of the kagome magnet is over 
70% of log 2. 

Most fundamentally, the single triangle is relatively much 
more degenerate than a tetrahedron, with 6 out of 8 (rather 
than out of 16) states being ground states. The triangle states 
have varying magnetizations of ±1 (whereas states obeying 
the ice rule have a unique magnetization, unless one tunes a 
field to a transition between magnetization plateaux^). 

To move towards the full triangle-based lattices, it is worth 
noting that for any Ising antiferromagnet, the 'hopping net- 
work' formed by neighbouring aligned spins can have a co- 
ordination of at most half that of z, the coordination of the 
underlying lattice - otherwise it would be energetically favor- 
able to flip the highly coordinated spin. The concept of the 
hopping network generalizes the loops on which the electrons 
hop on the planar and three dimensional pyrochlores. 

It now turns out that triangle and kagome lattices behave 
entirely differently from pyrochlore in both two and three di- 
mensions, on account of the nature of their frustrated ground 
state ensembles. 



1. Triangular lattice 

The ground-state degeneracy of the triangular lattice is im- 
mediately lifted by the addition of even a single hole. This 
result is entirely analogous to the frustrated Nagaoka theorem 
presented in Ref. [3jl in the context of the magnetic super- 
solid discussed there, and it is also connected to the triangular 
Bosonic supersolids22^£. 

As z = 6 for the triangular lattice, the hopping network of 
aligned spins no longer has coordination two as in the loops of 
the pyrochlore lattice. In fact, the coordination of a site of the 
hopping network no longer even needs to be uniform, so that 
there may be dangling or even isolated sites, as on Fig.[T3l(a). 
The coordination can range all the way from zero (for a spin 
surrounded by a hexagon of oppositely aligned spins) to max- 
imally z/2 = 3. 

The latter happens when the hole sits on a site experienc- 
ing zero net exchange field from its six neighbors. There ex- 
ists a unique state (pictured in Fig. [I3](b)) in which there is 
a network of three-fold coordinated sites. This state breaks 
translational symmetry by tripling the unit cell as well as time- 
reversal symmetry as it corresponds to a state with a magne- 
tization of a third of the saturated value. We have not stud- 
ied what happens to a finite doping but a Fermi liquid regime 
on the hexagonal backbone at low doping looks likely. This 
would imply a conducting state, with the possibility of ad- 
ditional low-energy excitations in the form of defects of the 
hexagonal backbone. 
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FIG. 12: Average number of filaments in 2D (left) and 3D (right), corresponding respectively to systems of sizes L = 10, 20, 40 and 
L — 6, 12, 18 (red squares, green dots, blue triangles) as a function of doping p. 



(a) (b) 




(c) (d) 




FIG. 13: Configurations on the triangular (top) and kagome (bot- 
tom) lattices respecting the antiferromagnetic frustrated constraints. 
(a) and (c) are random configurations, while (b) is the conducting 
hopping network minimizing the kinetic energy on the triangular lat- 
tice, and (d) is the insulating state maximizing the number of closed 
loops on kagome. Panel (d) also illustrates the arrangement of loops 
described in section |TV C 21 in the two dimensional kagome layers 
which are part of the three dimensional pyrochlore lattice. 



2. Kagome lattice 

The situation on the kagome lattice is different still. Its 
magnetic ground state ensemble is very short-range corre- 
lated, unlike that of the triangular lattice, which has algebraic 



correlations. With z = 4 for the kagome lattice, the hopping 
network can no longer branch but it can now have dangling 
links or isolated spins (coordination 1 and 0, respectively, as 
illustrated on Fig.[T3l(c)). 

The hopping networks thus consist of loop segments which 
need no longer close on themselves, minimally containing 
only one spin but not bounded above in the thermodynamic 
limit. The length distribution is, however, unlikely to contain 
long loops as this would only happen if each the magnetiza- 
tion of all the triangles along the loop segment has the same 
sign, at considerable cost in entropy. 

Energetically, it is of course again most advantageous to 
have closed loops, as electrons on them gain hopping energy 
— 2\t\. On the kagome lattice, such loops are readily con- 
structed. The shortest ones are obtained by arranging spins 
to be aligned around a hexagon. Their densest packing is ob- 
tained for a state which breaks translational symmetry, tripling 
the unit cell and incorporating a 1/3 magnetization as was the 
case of the triangular lattice above, see Fig.[T3](d). However, 
this state - which is the unique ground state at electron density 
p = 1/9 is now an insulating one - all hopping paths for the 
electrons close back on themselves after six steps. 
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Appendix A: Worm algorithm 

a. Worm algorithm without electrons. In our definition, 
a loop is made of nearest neighbour spins pointing in the 
same direction: it is uniquely defined and possesses an up 
or down flavor. On the other hand, a worm consists of al- 
ternating up/down/up/down/... spins; through each vertex of 
the premedial lattice (square in 2D or diamond in 3D); see 
figure[l4] The worm can randomly choose between two paths 
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energetically equivalent and eventually hit its initial position; 
the worm is then closed. Reversing all spins in the worm gives 
way to a new configuration in the Coulomb phase (figurefPfl). 
This method ensures both ergodicity and detailed balance to 
the algorithm in absence of electrons^. 

Worms, as defined above, are not self-avoiding, i.e. they 
can erase their own path. If a worm goes an odd number of 
times through the same site, the corresponding spin will be 
flipped during the Monte Carlo update; if it happens an even 
number of times, the spin shall not be flipped. 




FIG. 14: On the left, the same loop configuration as in figure [JJ 
The thick dashed green line represents a possible worm, as built by 
our algorithm; nipping the spins along this worm does not break the 
ice-rules (see on the right) and allows us to visit the configurational 
space of the Coulomb phase, ensuring ergodicity. 

b. Decorrelation. Since we are dealing with Ising spins, 
two random configurations roughly differ by one half of the 
system sites on average. Hence two configurations separated 
by k worm updates can be considered decorrelated if close to 
one half of the system sites are different: we arbitrarily chose 



45%. Previous works , has shown that a finite constant number 
of worm updates (~ 10) was enough to decorrelate the system 
in 3D, but an increasing number with system size L was nec- 
essary in 2D. This is related to Polya's theorem stating that 
a random walk (similar to the worm here) in 3D is transient 
(finite probability never to come back to the origin), while it 
is recurrent in 2D (it always comes back)2i 

c. Worm algorithm with electrons - entropic regime. In 
the entropic regime, i.e. at low doping, there is at most one 
electron per loop. A correct sampling can then be done by 
using the above worm algorithm free of electrons, and reject- 
ing all configurations with more electrons N e than loops Nf. 
However, the distribution of number of loops per configura- 
tion being Gaussian (see figure [3]), the density of electrons 
acts as a cutoff and above a certain threshold of the order of 
l/{£), almost all configurations are rejected and it becomes 
impossible to get good statistics. 

In order to take into account the influence of itinerant elec- 
trons in the entropy, we chose to weigh any loop configuration 
by the number of possible combinations to distribute N e elec- 
trons in Ng loops, namely the binomial coefficient C N e . 

d. Worm algorithm with electrons - energetic regime. In 
the energy regime, i.e. at intermediate or high doping, there is 
more than one electron per loop. The above method based on 
the worm algorithm without electrons become inefficient and 
finite temperature Monte Carlo simulations are necessary. 

Starting from a given configuration, a worm update is pro- 
posed. The loop histogram of the states before and after the 
proposed update are computed, and the corresponding elec- 
tron energies are calculated after filling the energy levels with 
N e electrons. The worm update is then accepted or not via a 
Metropolis argument. 
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